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Abstract 

In this work is obtained new approximation for the yield function Y for cosmic ray in- 
duced ionization in the Earth atmosphere on the basis of simulated data. The parameteri- 
zation is obtained using inverse nonlinear problem solution with afxy(analyze fx=y)-code. 
Short description of the methods is given. The found approximation is for primary proton nu- 
clei, inclined up to 70 degrees zenith angle. This permits to estimate the direct ionization by 
primary cosmic rays explicitly. The parameterization is applicable to the entire atmosphere, 
from ground level to upper atmosphere. Several implications of the found parameterization 
are discussed. 
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1 Introduction 



Presently it is known that the Earth is hit by elementary particles and atomic nuclei of very large 
energies and in wide energy range, this is the cosmic ray radiation. The fluxes variate, from 10 4 
m~ 2 s^ 1 at energies 10 9 eV to 10~ 2 kmT 1 yr~ l at energies 10 20 eV. The cosmic ray intensity 
is approximatively expressed with (1), where E is the total particle energy per nucleon in GeV 
and a = —2.1 is the spectral index. The majority of these particles are protons. The primary 
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cosmic rays penetrate upper atmosphere, secondary cosmic rays are produced by the interactions 
of primary rays in the atmosphere. 

I n (E) a \.%{-^-^) a nucleons. cmT 1 .s~ l .sr~ l (1) 

The abundances are approximately independent of energy, at least over the dominant energy 
range of 10 MeV/nucleon through several GeV/nucleon. By mass about 79 percent of nucleons 
in cosmic rays are free protons, and about 80 percent of the remaining nucleons are bound in 
helium. 

There are three different types of cosmic rays: galactic cosmic rays, solar cosmic rays and 
anomalous cosmic rays. Most galactic cosmic rays are accelerated in the shock waves of su- 
pernova remnants. Because of their deflection by magnetic fields, galactic cosmic rays follow 
convoluted paths and arrive at the top of the Earth's atmosphere in practice uniformly from all di- 
rections. Galactic cosmic rays are the most typical cosmic rays, and their flux in the solar system 
is modulated by the solar activity: enhanced solar wind shields the system from these particles 

Q. 

We know that the galactic cosmic rays create the ionization in the stratosphere and troposphere 
and also in the independent ionosphere layer at altitudes 50-80 km in the D region 0. This 
ionization is a result from the impact of the secondary cosmic ray electromagnetic, muon and 
hadronic components on the planetary atmosphere. The cosmic ray induced ionization is an 
important factor of space influences on atmospheric properties. 

In general the variations of cosmic ray induced ionization are caused by solar activity varia- 
tions, which modulate the cosmic ray flux in interplanetary space, and changes of the geomagnetic 
field, which affects the cosmic ray penetration in the atmosphere and their access to Earth. The 
changing solar activity is responsible for the variation of solar wind, respectively cosmic rays. 
The solar wind reduces the flux of cosmic ray reaching the Earth, since a larger amount of en- 
ergy is lost as they propagate up the solar wind. Since cosmic rays dominate the ionization, an 
increased solar activity will translate into a reduced ionization. 

On the basis of satellite data analysis was shown that cloud cover varies with the variable 
cosmic ray flux reaching the Earth 01 31. Over the relevant time scale, the largest variations arise 
from the 11 -yr solar cycle, and indeed, this cloud cover seemed to follow the cycle and a half 
of cosmic ray flux modulation. Specifically was shown that the correlation is primarily with low 
altitude cloud cover jH. 

In this connection a detailed model of the cosmic ray induced ionization will be a good basis 
for a quantitative study of different mechanisms affecting Earth's atmosphere. To estimate the 
cosmic ray induced ionization it is possible to use a model based on an analytical approximation 
of the atmospheric cascade flU or on a Monte Carlo simulation of the atmospheric cascade |[71 
[H. A key issue, which allows to estimate the cosmic ray induced ionization for given location, 
altitude and spectrum of cosmic rays is the use of ionization yield function Y (2) which is defined 
according [[91 

Y(x,E) = AE(x,E)^--!-Q. (2) 

where AE is the deposited energy in layer Ax in the atmosphere and Q. is a geometry factor, 
integration over the solid angle with zenith of 70 degrees. Afterwards the ion pair production q 
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by cosmic rays following steep spectrum is easy calculated according the expression: 

q(h,X m ) = J D(E,X m )Y(h,E)p(h)dE (3) 

where D(E,X m ) is the differential primary cosmic ray spectrum at given geomagnetic latitude, 
Y is the yield function according (2), p(h) is the atmospheric density in [g.cm~ 3 ]. 

The ionization yield function Y depends only of assumed physical models for cascade pro- 
cesses in the atmosphere, the atmospheric model and the type of the primary particle. Therefore 
having a convenient parameterization for ionization yield function Y and parameterized differen- 
tial energy spectrum of galactic cosmic rays at the Earth's orbit it is possible to estimate cosmic 
ray induced ionization in different locations and conditions. In this work we use the previously 
obtained ionization yield function Y [fTOlfTTfl for primary protons. 

The major contribution of cosmic ray induced ionization in the atmosphere is due to particles 
having energy till 1 TeV, taking into account the steep spectrum of primary cosmic rays (1). This 
is the reason to deal in this paper with primary protons with maximal energy of 1 TeV. 



2 How to investigate a given nonlinear system 

Many problems in physics, applied mathematics even and pure mathematics lead to solution of 
nonlinear system of equation. 

To analyze nonlinear equation, we write it in the form 

fx = y, f:D f cR n ^M."\xeR'\yem m . (4) 

To analyze iteratively nonlinear systems involves two related classes of problems [fT2l - heuris- 
tic investigation of nonlinear systems issuing from not-refined mathematical models and auto- 
matic solution of streams of one-type nonlinear systems of equations. 



2.1 Main iteration procedure in a/xy-code[12J 

A powerfull tool for an heuristic analysis of systems of nonlinear equations (4) is q/ry-code 
(analyze fx=y).In the case n < m as well in non-trivial case «>mwe employ regularized Gauss- 
Newton type iterator 

x £ D f , x° £ Df, £° > 0, (f' T (x k )f(x k ) + (e k + a k )I){x k+l -x k ) = (5) 

-f fT (x k )(fx-y)-a k x),k = 0,l,... . 

The regularizators £ k (of the process) and a k (of the problem) will be defined below. 

In the right hand of (5) original problem (4) is transformed by Gauss manner Fx : = / T (x) (fx — 
y) and additively regularized by von Neumann-Tichonov manner. The space W we consider as a 
two different normed elementwise equal spaces R^, and RJ}- 

The space is normed by norm ||a||<x> = max|a,|, (a = [a\,a2, ...,a„] r £ IR' 1 ) and 

space is normed respectively by L2 norm ||c||2 = (E/Li a f) 1 - 
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The corresponding matrix spaces [R£,,R£J and [IR^IR^] are normed respectively by [jAjjoo = 
max L'/=i \ a ij\i A = {«(/} G [R^R^,] and by ||A||2 = hn, where X m is the maximal spectral num- 

ber of the matrix A [Tf3l . Note, the norm ||A||oc is well defined. It is simultaneously completed 
and consisted to the vector norm ||a||«>. 

The iterator (5) produces two different approximating sequences subjects to the used regular- 
izators. 

To solve an ill-conditioned [14] problem (4) we use the process (5) in the R£, space with 
autoregularizator of the process [fT51 [T6l 

e k =\{j{x k Y + AN P k -x k ^, (6) 

where x k = ||/ /7 V)/V)IU p k = \\f' T (x k )(fx k -y)^ and N = const > 0. The behavior of 
the autoregularizator is accordant semilocal convergence lfT6l of the process (5). As a rule the 
convergence domain expands with increasing the initial value of the regularizator e . 

To solve a singular problem [141(4), we use the process (5) in the space Rj w i m autoregular- 
izator (6) with 

T* = A* (7) 
is the smallest non-negative eigenvalue of the matrix f' T \x k )f'(x k ) and 

p k =\\f' T (x k )(fx k -y)-a k x)\\ 2 , a k = rexp(-a(k-l)), y,o = const>Q. (8) 

Presently in the case of singular problem we have not good semilocal theory for the conver- 
gence of the process (5)-(8). In the case of when convergence persists to a point x*, the value 
||x* — x|| 2 is minimal. Usually one takes x = 0. However the proposed iterator is fully applicable 
and was used for analysis of large diversity of problems. 

The described main iteration procedure starting from Levenberg-Marquardt years till present 
days was applied extremely successful. 



2.2 From local root extractors to global iterative processes 

In order to find all solutions of Eq. (HJ) in the domain Df the vector Fx 1 is repeatedly multiplied 
by the local root extractor 



e j (xJ J ^ = - 



l-exp(-||x-JcO')|||)' 

in which is the y'-th solution of Eq. ©. In the repeated solutions of the transformed problem 

F J x:= (jJe j (x,x U) )jFx = 0, J>1, (9) 

the process © is executed with a new Fx := F J x. For every solution the process © is started 
many times with different x° and e . Each time when J increases the derivatives f'(x ) are 
computed analytically and the matrices /' (x k )f'(x k ) are adaptively scaled. 



4 



To analyze globally nonlinear problem (4) means to find all solutions with estimation of them- 
selves isolation in the domain Df. For this purpose was build the mentioned above local root 
extractor. 



3 The parameterization of ionization yield function Y 

The ionization yield function Y gives the number of ion pairs, produced in 1 g of the ambient 
air at a given atmospheric depth by one particle of the primary cosmic ray radiation with given 
kinetic energy per nucleon. 

The approximation is obtained following procedures similar to described in lfT71[T8l . In this 
case we take the logarithm of the problem, which permits to approximate large diversity of distri- 
butions with similar shapes, having different amplitudes with one model function with different 
parameters [fT9ll20ll2TTl. 



3.1 Approximation of ionization yield function Y 

Usually the proposed approximation are carried out with large diversity and classes of functions, 
some very complicated. It exist powerful class of approximations - fractional-rational functions 
applied in our case. As a result we obtain description of the simulated data and fast convergence 
of the iteration process. For example we utilize simple fractional-rational function 

where h is the atmospheric depth and a,b,c,d,e and f are parameters. The fits for ionization 
yield function Y are presented in Fig. 1-8 for protons with energies 500 MeV, 1 GeV, 5 GeV, 10 
GeV, 50 GeV, 100 GeV, 500 GeV and 1 TeV In the figures with solid black squares are shown 
the obtained approximations and with open circles the simulated data. During the solution of 
the inverse problem we take logarithm of problem and scaling of the argument. This permits to 
obtain parameterization with small uncertainties and fast convergence of the process. In this case 
we solve ill-conditioned problem with large number of condition 10 17 . The solution is obtained 
after only 100 iterations steps. In this case the normalized x 2 varies between 0.13 for 500 MeV 
protons and 0.001 for 1 TeV protons. 
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Figure 1 : Approximation of ionization yield function Y for 500 MeV primary protons 




Figure 2: Approximation of ionization yield function Y for 1 GeV primary protons 
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Figure 3: Approximation of ionization yield function Y for 5 GeV primary protons 
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Figure 4: Approximation of ionization yield function Y for 10 GeV primary protons 
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Figure 5: Approximation of ionization yield function Y for 50 GeV primary protons 
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Figure 6: Approximation of ionization yield function Y for 100 GeV primary protons 
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Figure 7: Approximation of ionization yield function Y for 500 GeV primary protons 
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Figure 8: Approximation of ionization yield function Y for 1 TeV primary protons 
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3.2 Main results and parameters 

The quality of the parameterization is the same in all cases. 

The general aim is in one hand the precise description of the position of the Pfotzer maximum 
and on the other hand good description in the low atmosphere. 

In the case of ionization yield function for 500 MeV protons Fig.l the end points of the distri- 
bution are equally weighted. Our approximation is chosen to take the mean values at distribution 
tail. 

In the remain cases the main difficulties are connected with precise description of the Pfotzer 
maximum. As example for the ionization yield function for 1 GeV and 5 GeV protons on observes 
very good coincidence between approximation and simulated data (see Fig. 2 and Fig.3). The 
approximated ionization yield function for 10 GeV protons coincides with simulated data (Fig .4). 

In the case for 50 and 100 GeV, the parameterization gives slight increase of ion pairs in 
Pfotzer maximum (Fig. 5 and Fig.6). 

Finally in the case of ionization yield function for 500 GeV and 1 TeV protons the approxi- 
mation coincides once more with simulated data (Fig. 7 and Fig. 8). 

In all cases the approximation in practice coincides with simulated data in all other regions 
of the atmosphere,especially in the lower atmosphere. The parameters of the approximation are 
presented in Table. 1. 



Table 1 : The parameters of the approximation 



Energy 


a 


b 


c 


d 


e 


f 


500 MeV 


-6.09548 


3.14499 


0.31817 


2.76023 


0.37625 


0.06494 


1 GeV 


-2.51931 


3.81687 


0.78784 


1.0136 


0.64061 


0.15493 


5 GeV 


-4.55444 


6.69186 


0.47114 


0.2117 


1.0687 


0.09188 


10 GeV 


-0.89284 


4.01468 


0.19579 


0.30881 


0.5975 


0.03812 


50 GeV 


1.82327 


2.46817 


0.15813 


0.55512 


0.31109 


0.02995 


100 GeV 


2.28484 


1.99474 


0.13683 


0.56878 


0.22815 


0.0255 


500 GeV 


0.99557 


2.65916 


0.10061 


0.29285 


0.30755 


0.01891 


1 TeV 


0.03717 


2.72562 


0.24925 


0.13949 


0.27453 


0.04701 



The approximation parameters are easy to fit with exponential or power low functions. This 
permits using the expression for ionization (3) to estimate the cosmic ray induced ionization at 
given location and altitude. 

In addition, it is important that the approximation (10) is analytically integrable , which gives 
the possibility to estimate the total atmospheric ionization explicitly. 

4 Discussion 

The obtained parameterization of ionization yield function Y permits easily to compute the ion- 
ization rates from cosmic rays in the Earth atmosphere for given location and conditions, instead 
of using counting rates from different devices as a proxy. Moreover the parameterization allows 
to estimate variations of the cosmic ray induced ionization caused by the variable solar activity 
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and as result to evaluate several atmospheric effects of cosmic rays, on different timescales and 
under different helio spheric conditions. 

The cosmic rays may affect climate, and are probably important climate driver. Recently, it 
was shown in Il22ll that the variations in the amount of cloud cover follow the expectations from 
a cosmic -ray cloud cover link, specifically for low altitudes [5]. It was shown that the relative 
change in the low altitude cloud cover is proportional to the relative change in the solar-cycle 
induced atmospheric ionization at the given geomagnetic latitudes. Namely, at higher latitudes 
the ionization variations are about twice as large as those of low latitudes. 

Above 100 percent saturation, the phase of water is liquid and it will not be able to condense 
unless it has a surface to. Therefore for formation of cloud droplets the air must have cloud 
condensation nuclei. Changing the density of these particles, the properties of the clouds can 
be varied. Having more cloud condensation nuclei, the cloud droplets are more numerous but 
smaller, this tends to make whiter and longer living clouds. The suggested hypothesis, is that in 
regions devoid of dust, the formation of cloud condensation nuclei takes place from the growth 
of small aerosol clusters, and that the formation of the latter is governed by the availability of 
charge, such that charged aerosol clusters are more stable and can grow while neutral clusters can 
more easily break apart. 

If this process is dominant, charge and therefore cosmic ray ionization would play an impor- 
tant role in the formation of cloud condensation nuclei. With this in mind assuming our param- 
eterization (10) for ionization yield function Y and relation (3), as well the parameterization for 
differential cosmic ray spectrum [23J we estimate the ionization rates for several cut-off rigidities 
(Fig.9). 
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Figure 9: Ion pairs for Proton induced showers for 1, 5, 9 and 15 GV rigidities 
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The presented results demonstrate realistic cosmic ray induced ionization profile using av- 
erage primary cosmic ray spectrum during moderate solar activity. As it is seen these results 
are applicable for different geomagnetic latitudes, respectively different locations. The obtained 
results are in good agreement with the profiles described in [|9l[TTI. 

It exists the hypothesis that the intrinsic variation in the cosmic ray flux are clearly evident in 
the geological paleoclimate data. Within the determinations of the period and phase of the spiral- 
arm climate connection, the astronomical determinations of the relative velocity agree with the 
geological sedimentation record for when Earth was in a hothouse or icehouse conditions lF24l . 
Moreover, it was found that the cosmic ray flux can be independently reconstructed using the so 
called "exposure ages" of Iron meteorites ll25l . 

Thus it is very important, in one hand to investigate the hypothetic influence of cosmic rays on 
the Earth's climate trough the mechanism cosmic rays induced ionization-condensation nuclei- 
cloud condensation nuclei-cloud cover, and on the other hand to study the impact of intrinsic 
variation in the cosmic ray flux during the motion of our planet in the Galaxy in attempt to 
estimate and clarify their contribution. 

5 Conclusions 

In this work is presented new convenient parameterization of the yield ionization function Y of 
primary cosmic rays into the Earth atmosphere. The parameterization was obtained on the basis 
of inverse problem solution using afxy-code. 

The importance of cosmic ray induced ionization is widely discussed, precisely the mecha- 
nism related to clod cover formation. 

Our results are important for precise estimation of the ionization profiles in the atmosphere, 
when one deal with proton nuclei from primary cosmic ray. 
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